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Abstract: We address some recent concerns about the absence of the Boulware-Deser ghost 



in the Stiickelberg formulation of nonlinear massive gravity. First we provide general argu- 
ments for why any ghost analysis in the Stiickelberg formulation has to agree with existing 
■j . consistency proofs that have been carried out without using Stiickelberg fields. We then 

J> , demonstrate the absence of the ghost at the completely nonlinear level in the Stiickelberg 

en 

00 

fN| | field and the associated secondary constraint that eliminates its conjugate momentum are 

computed explicitly, confirming the consistency of the theory in the Stiickelberg formulation. 
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1 Introduction and statement of the problem 

General relativity describes a massless spin-2 field g^ u . Constructing a mass term for g uv 
which is reparametrization invariant requires introducing another "metric", say fu V - Then 
a potential can be constructed as a function of g^fxv, and one can write a massive gravity 
action in the form, 



S = M 2 V J d^x^-g [R(g) - 2m 2 V(g~ l f)] . (1.1) 



Alternatively, this describes a massive spin-2 field in fixed background metric fau- 
lt has long been known that for a generic V such a theory is inconsistent, containing 
the Boulware-Deser ghost [1, 2]. Only recently potentially consistent massive gravity theories 
were proposed mostly based on a perturbative analysis in the Stiickelberg formulation of 
the theory [3, 4]. These were then generalized and shown to be ghost free in a completely 
nonlinear analysis without using the Stiickelberg formalism [5-8]. In spite of these proofs, 
concerns about the consistency of massive gravity in the Stiickelberg formulation has arisen, 
both at the perturbative [9, 10] as well as nonlinear levels [11, 12]. At the perturbative level 
the problem has been addressed in [13-15]. Here these issues will be addressed and resolved 
at the completely nonlinear level. 

We start with stating the problem. The simplest choice for f^ u in (1.1) is the flat metric, 

, d<j> A 34> B 

^ = ^AB-^. (1.2) 

The four scalar fields ensure reparametrization invariance of V(g _1 f) with flat / and 
can be gauged away by a reparametrization x = <f> (x) to set f^ v = rj^. These are the 
Stiickelberg fields that inherit their dynamics from the form of V. It was pointed out in 
[16, 17] that the consistency of the <j> theory is correlated to the existence of the Boulware- 
Deser ghost. This allowed for a simpler analysis of the ghost which was otherwise an involved 
problem. In particular, one could work in a "decoupling limit" that ignores the nonlinear 
dynamics of g^ u and retains only the <fi . By the argument of [16], this retains also a part of 
the ghost information that survives the decoupling limit. The potentially ghost-free poten- 
tial V{g~ l f) based on (1.2) was first constructed [3, 4] using a perturbative analysis of the 
Stiickelberg fields in the decoupling limit. It was established to be ghost free to linear order 
in hu V = g^ u — rjuu . A fourth order analysis was also carried out [4] . 

The subsequent development in the field did not rely on the Stiickelberg formulation. 
The nonlinear demonstration of the absence of ghost in the flat / theory of [4] was carried 
out in the unitary gauge / = r\ [6]. Subsequently, these actions were extended, first, to any 
arbitrary non-dynamical fn, v [7, 8], and then to a dynamical f^ u [8, 18] and shown to be ghost 
free at the complete nonlinear level, without a need for gauge fixing. This established that 
such theories contained the right constraints to eliminate a propagating ghost. 

On the other hand, while the Stiickelberg analysis, based on studying the dynamics of the 
<f> (1.2), proved very powerful in the decoupling limit, extending it beyond this limit led to 



some speculation that the ghost may return at higher orders in h^ u [9, 10]. This concern was 
addressed in [13] where it was asserted that the constraints that eliminated the ghost [6-8, 18] 
must also arise in the Stiickelberg formulation of the same theory. This was demonstrated 
perturbatively at lowest orders in h^ u (as well as exactly in 2 dimensions) and argued to be 
extendable to higher orders. Recently, a nonlinear analysis of constraints in the Stiickelberg 
setup has been performed [11, 12] and seems to indicate that the required constraints may 
not exist in the Stiickelberg formulation. 1 If true, this would be in contradiction with the 
proofs of [6-8, 18] as well as the arguments in [13]. 

Due to these lingering doubts about the consistency of massive gravity in the Stiickelberg 
formulations and in order to conclusively resolve this confusion, we reconsider this set up. 
Our main results are summarized below. 

• With simple arguments based on general covariance, we show that the Stiickelberg 
formulation of massive gravity must agree with the existing consistency proofs already 
obtained in the standard formulation. 

• We perform a Hamiltonian (ADM) analysis of the minimal model of massive gravity in 
the Stiickelberg formulation and explicitly obtain 2 constraints that remove both the 
ghost and its canonical momentum. The analysis is performed at the fully nonlinear 
level and for arbitrary metric f^ u . Our results extend those of [13] which, for a flat 
f^ u , found the corresponding constraints exactly in 2 dimensions and perturbatively in 
4 dimensions. This explicit nonlinear proof should settle the issue of absence of ghost 
in the Stiickelberg formulation of massive gravity. 

In section 2, we discuss the consistency between the Stiickelberg and standard non-Stiickelberg 
formulations of massive gravity. In section 3, we perform the consistency analysis directly in 
the Stiickelberg formalism and prove the absence of ghost at the nonlinear level. 

2 Consistency of the Stiickelberg and non-Stiickelberg formulations 

The absence of ghost in massive gravity has been proven in the standard, non-Stiickelberg, 
formulation of the theory [6-8]. It is easy to show, on general grounds, that this also implies 
the absence of ghost in the Stiickelberg formulation. 

In the standard formulation, the fu V in (1.1) is a non-dynamical tensor that ensures 
general covariance, while the equations of motion are obtained by varying the action with 
respect to g^ alone, 

R^u ~ \g^R - 2m 2 V^ u = , (2.1) 

where, V^ u = j-^j — -^g^V . The Bianchi identity for the Einstein tensor then implies, 

vnv = °- (2.2) 



x This is different from the objection raised in [19] and addressed in [8]. 



On the other hand, in the Stiickelberg formulation of (1.1), one writes [16], 

d(j) A - d(p B 

Jab are given scalar functions of x 1 ^, but the Stiickelberg fields <\> are now treated as dy- 
namical. Hence along with the g^ equation above, one also obtains the <fi equations of 
motion, 

£-■ 

However, general covariance implies that (2.4) is already contained in the g^ v equations (2.2). 
To see this, consider infinitesimal reparametrizations 5x^ = £^ under which, 5g^ u = — 2V^£^ 
and 5<p = —^d^cj). The invariance of the action implies, 



8S = -Ml J d 4 x^g- (R^ - \ 9ilv R - 2m 2 V^) SgT + 2m 



2 SV xaA 



0. (2.5) 



For the above variations, integrating by parts and using the Bianchi identity then gives, 

v/ %-=£z^- ( 2 - 6 ) 

Since <f> are non-singular coordinate transformations of x^, d^cfr is an invertible matrix. 
Then (2.6) implies that (2.4) is equivalent to (2.2), as can be explicitly checked [5]. This 
is a consequence of general covariance and shows that the Stiickelberg fields do not lead to 
extra equations beyond the metric equations of motion. Adding matter does not change this 
argument. 

Now, the proof of absence of ghost in massive gravity in [6-8] only involves the g^ u 
equations of motion. In particular, it does not involve the f^ v equations nor gauge fixing 
(except for the unitary gauge analysis in [6]). Therefore, in the final action for the five 
physical components of the massive spin-2 field, obtained on imposing the constraints and 
eliminating the ghost, one may express / as in (2.3) and obtain the <f> equations. But as 
discussed above, these are already contained in the g equations, showing that the Stiickelberg 
formalism cannot be inconsistent with the ghost analysis of [6-8]. 

After this general argument, let us briefly review the ghost analysis for the simplest of 
the massive gravity actions, the minimal model of [5], following [7, 8], 



S = M 2 V J d A x^~g [R(g) - 2m 2 Tr( Vg^f)] . (2.7) 



Here, \J g 1 f is a square-root matrix defined such that \JE\JE = E. To obtain the physical 
degrees of freedom one uses the ADM parametrization, 

_ / N 2 + N i lijNJ N \ J_/-l N* \ 

3^~\ N . 7 . I > 9 - N2 \ NJ N 2 fy i j _ N i N j ■ W 



and writes the action in the Hamiltonian formulation, 

S = Ml j d 4 x frijij + iV.Ro + N l Ri - 2m 2 Ni/d^ V min (N, N\j, /)] . (2.9) 



where, V m in stands for Tr(y/g~ 1 f) expressed in the ADM parameterization and, 

Ro = Vl ?R + {W ~ Tf'-Kij) lyfi] , Ri = 1y/l Vj {^/^f) . (2.10) 

In the Hamiltonian formulation, the action contains 6 potentially propagating modes jij and 
their canonically conjugate momenta -n 13 = 5S/5jij. Of these, 5 conjugate pairs describe the 
massive spin-2 graviton, while the sixth one is the Boulware-Deser ghost. The N and Ni 
have no canonical momenta and are non-propagating. All these fields give rise to equations 
of motion. The crucial point is to show that these equations contain the two constraints that 
can eliminate the ghost and its conjugate momentum. There are no equations of motion for 
the spectator fields f^. 

Of the four N and N l equations of motion, three combinations determine the N l in terms 
of N, jij and tt ij (the explicit solutions are obtained for three functions n l (N,N l ) in terms 
of jij and 7r JJ ) . The fourth combination becomes the Hamiltonian constraint on jij and ir lJ , 

C(j ij ,^) = 0, (2.11) 

where C does not involve time-derivatives of its arguments. The preservation of this constraint 
by the time evolution of the system requires dC/dt = 0. Eliminating all time derivatives of 
fields using their dynamical equations leads to a second constraint, 

C 2 (jij,^) = 0. (2.12) 

These two constraints eliminate the ghost field and its conjugate momentum. Finally N is 
eliminated by the equation following from the preservation of C 2 in time. All these equations 
are part of the g^ v equations of motion and leave behind a theory for the 5 physical modes 
of a massive spin-2 field. 

The final theory also contains L which was treated as a spectator all along. If this is 
now expressed in terms of the Stiickelberg fields 4> (2.3), we know that their equations of 
motion are already part of the g^ u equations and there is no inconsistency. In addition, general 
covariance should be used to eliminate four gauge modes. Below, we perform a Hamiltonian 
analysis directly in the Stiickelberg formulation of the theory. 

3 ADM analysis in the Stiickelberg formulation 

In this section we will perform a Hamiltonian analysis of massive gravity in the Stiickelberg 
formulation and derive the constraints needed to eliminate the ghost. This explicitly answers 
the concerns raised in [9-12], about the consistency of massive gravity in the Stiickelberg 
formulation and extends the perturbative arguments of [13] to the completely nonlinear theory. 



3.1 Hamiltonian form of the action with Stiickelberg fields 

Here we explicitly focus on the minimal massive gravity action (2.7) in the Stiickelberg for- 
mulation with fyw given by (2.3) and reconsider the analysis of this model in [11, 12]. To 
avoid working with the square-root matrix, we follow [20, 21] to recast the action as, 

S = Ml,Jd i x^[R{g)-m 2 [^ A A + ^- l ) A B A B A ]] , (3.1) 

where, 

A B A = d^ B g^d^ C fcA. (3.2) 

Note that since A = A bj is symmetric, only the symmetric piece of <&ab — Iac^b 
appears in the action. We therefore treat §ab as symmetric in the following. Solving the 
$ B equation of motion and plugging the solution back into the action gives back (2.7). 

To rewrite the action (3.1) in the Hamiltonian formulation, we summarize the analysis of 
[11, 12] until our conclusions diverge, and mainly adapt their notation to facilitate comparison. 
In the ADM parametrization (2.8), the matrix A B A of (3.2) can be written as, 

A B A = -Vn^Vn^ JcA + V B A , (3.3) 

where, 



V n A = -(d cf> A - N l d i( p A ) , V B A = -yVdi^drfCfcA • (3.4) 



1 

N 
In this notation, the canonical momentum conjugate to eft is, 



S = Ml I d 4 x 



p B = 2m 2 Ml^i{<S>- l ) AB V n (k A . (3.5) 

Hence, the action in the Hamiltonian formulation becomes, 2 

^j^+pA^-n] , (3.6) 

with the Hamiltonian density given by, 

n = -NR - WRi + NH s f + N l Uf . (3.7) 

The first two terms are familiar from general relativity, whereas the rest read as, 

«5* = V^mW(* a a + (<s>-y B v B A) + ± m2 l M 2^ * ab vavb , 



lA 



Uf=p A di4> A . (3.8) 

Before we set out to determine the physical content of the theory, let us begin by listing 
its field content in the Hamiltonian description: 



2 To avoid lengthy expressions, we mostly work with the action principle in the Hamiltonian formulation. 
Only in section 3 we work with Poisson brackets in a minimal way. 



fields 


components 


canonical momenta 


N 


1 


- 


Nt 


3 


- 


Hj 


6 


7T« 


cP A 


4 


PA 


$ A B 


10 


- 



Thus there are 10 potentially propagating modes given by 7y and the Stiickelberg fields </> , 
while 14 further fields are non-propagating. There are also 4 gauge invariances. 

A naive counting of the field components and non-dynamical equations of motion (ones 
without time derivatives) derived from (3.6) may give the impression that the theory contains 
the Boulware-Deser ghost: At first sight, the 10 equations of motion for the $ b depend on 
$ B an d will therefore determine this matrix rather than serve as constraints on other fields. 
Since the action is linear in iV and N l , their equations of motion turn into 4 non-dynamical 
equations for the remaining variables. But unlike general relativity, now these also contain 
the pa and <fi and can be solved, for example, for the pa (after gauge fixing the (f> ) rather 
than impose a constraint on the Boulware-Deser ghost contained in 7j,-. Hence, it seems that 
the theory may not have the required constraints to reduce the number of propagating modes 
of jij below 6, as concluded in [11, 12]. 

However, as we will demonstrate in the following, the 10 equations of motion for 3> b 
depend only on 9 independent combinations of the matrix elements $ b- One combination 
remains undetermined and the corresponding equation of motion will give an additional con- 
straint on the remaining variables instead, which will remove the ghost. Below we obtain this 
and the associated secondary constraint. 

3.2 The first constraint 

Varying the action (3.6) with respect to Q AB gives the equations of motion [12], 

1 



^AB = VlM£m 2 (f AB ~ (^- 1 )acV ou {^- 1 ) db ) + 



4m 2 M|^/7 



PaPb = 0. 



(3.9) 



To arrive at these equations we have divided by N which is non-zero since g^ is invertible. 
Since ^ ab is symmetric, the naive expectation would be that (3.9) provides 10 conditions 
that determine all 10 components of the symmetric ($ _1 ) . 

However, the crucial observation is that the 4x4 matrix V AB = j'^dicj) dj(j) B always has 
rank 3 since it is composed of 3 x 4 and 3x3 matrices each of maximum rank 3. This also 
implies that W = &~ 1 V&~ 1 has rank 3 and hence cannot depend on more than 9 independent 
combination of the 10 &ab- To see this, note that, being a symmetric rank-3 matrix, W can 
be diagonalized by an appropriate orthogonal transformation O to, 



W B = O t WO = diag{wi,w 2 ,w 3 ,0} 



(3.10) 



where W-q always has one zero eigenvalue. Since the orthogonal matrix O has six independent 
parameters, on inverting the transformation we see that W depends on at most 9 indepen- 
dent parameters. This shows that the 10 equations (3.9) depend only on the 9 independent 
combinations of <& AB that appear in $ _1 y$ _1 and can only determine these combinations. 
One combination of the equations therefore cannot fix $ , rather it constrains 7y and pa- 
The constraint hidden in (3.9) is extracted by multiplying on both sides with $ to get, 

As noted above, V is a rank-3 matrix and so the left-hand side also has to have rank 3. 
By definition, $ is invertible and thus has rank 4. Therefore the matrix 

CAB^fAB + J^^—^PAPB (3.12) 

is constrained to have rank 3 by equation (3.11). This implies the constraint C = det C = 0. 
Since paPb is a rank-1 matrix, det C is particularly simple and the constraint becomes, 

c " / ( /AB $r +1 ) =0, (3,13) 

where / = det(fAB) and a = 2m 2 Mp. This is a single constraint on 7™ and pa and can be 
used to determine 7 = det (7^) in terms of the pa- Together with the equations obtained 
from varying the action with respect to N and N l (which can be solved for the pa on gauge 
fixing the <fi fields), this gives a single constraint on the 7^ and the 7r lJ that can remove the 
ghost field. The associated secondary constraint will be obtained in the next subsection. 

To compare with existing results, in [13] the analogue of (3.13) was obtained in 2 dimen- 
sions and also perturbatively in 4 dimensions to linear order in h^ v = g^u — V/j-f (i n this case, 
for a different choice of the potential). Here we have derived the constraint in the minimal 
model of massive gravity in Stiickelberg formulation at the fully nonlinear level. The proof 
can be generalized to any dimension and, unlike previous studies, is valid for general Jab- 

We now obtain a result that will be used in the next subsection. Let us write (3.9) in the 
form $ i&ABP = 0. On imposing C = and using V given in (3.4) this becomes, 

V CD (<S>- l ) D A p A = d t ct> c ^d j( t> D {<S>- l ) D A pA « 0. (3.14) 

In words, ($~ 1 )a Pc is th e eigenvector with zero eigenvalue of the rank-3 matrix Vb ■ These 
equations imply a result that will be used later, 

7 fei a i D ($- 1 )z)V«O. (3.15) 

The last step can be justified explicitly. The map di<p can produce any 3- vector Bi as a map 
from some 4- vector be- In particular, one can always find three 4- vectors b^ (a = 1,2,3) that 
map to 3 linearly independent 3- vectors Bf, 

B? = di(j) c 'b a c , a = 1,2,3. (3.16) 



Then, Bi a regarded as 3 x 3 matrix is invertible (easiest to verify in a unitary gauge <jP = x c ). 
Now multiplying (3.14) with {B~ v ) k a {b a ) c gives (3.15). 

3.3 The secondary constraint 

In order to eliminate both the ghost and its conjugate momentum, a second constraint is 
needed. This should arise as a secondary constraint that preserves (3.13) under time evolution, 
i.e., £2(7,71") oc dC/dt pa 0. It is most efficient to compute dC/dt using Poisson brackets, 3 

§ = {C , Ft*} • (3-17) 

Below, we compute this and verify that it indeed leads to a constraint. 

An economical way of going over from the action (3.6) to the Poisson bracket formulation, 
is to regard only (7^ , vr*- 7 ' ) and (0 ,pa) a s canonical pairs and obtain all non-dynamical 
equations in terms of Lagrange multipliers. This requires introducing the Hamiltonian, 

H tot = f d 3 x [N(H% -R + r AB V BA ) + N\Uf - Ri)] . (3.18) 

Note that as compared to (3.7) this contains the extra Lagrange multipliers r to obtain 
the <& AB equations of motion N^ba (3.9). Of course, the theory can also be extended by 
introducing extra momenta and Lagrange multipliers as in [12]. But that will not affect the 
computation here. The r are determined by the new <&ab equations of motion. 
To compute the Poisson brackets we need the following identities, 

{p A {y),d xl (j) B {x)} = -8 A d yi 8(x - y) , {pa(v),Pb(x)} = 0, 
ifBc(y),PA(x)} = 6 -j^-, {THv),^*)} = 0, (3.19) 

together with (where V is the 7^ compatible covariant derivative), 

5RAx) 

5 Jki z \ = -(7ifc(z)V 2J + -fij(x)V zk )5(x - z) , 

Hi! = -^( 7 i*(*M*) - 2, M X))S(X - z) . (3.20) 

For {C, -Htot} to give a constraint on 7 and ir, it should not depend on N,N' 1 and the single 
component of <& AB not determined by (3.9). Using the above identities we find, 

{C{y),Hf{x)} « -2f(y)V yl 5(x-y) « {C(y), R t (x)} ■ (3.21) 

Hence the term proportional to the N l in H to t does not contribute to C2 on the constraint 
surface. The remaining terms are, 

{C(y), H tot } = J d 3 x N(x){C(y) , U^{x) + T A B (x)* B A (x) - R (x) } . (3.22) 



3 The Poisson bracket is denned as {f(x),g(y)} = / d 3 z ( //igj //^ - g) £$ ) , where q a and Po 
stand for the variables <f> , jij and their canonical momenta pa, tt 1 -' , respectively. 



Each bracket can be evaluated on the constraint surface using (3.19) along with, 

{ PA (y), V c D (x)} = - 1 i \x)d xi( f> B (x) [{5 c A f BD + 5 c B f AD )d y] + d xJ c5 -0-] J(x - y) , 

{C(y),R (x)} « fn-y-^Six - y) . (3.23) 

The coefficient of the term proportional to djN vanishes by virtue of (3.15). It is then 
straightforward to see that the result is proportional to N(y) which therefore appears as an 
overall factor in {C,H tot } and can be divided out. Then the secondary constraint reads, 

C 2 = N-^CiylHu*} « ~T 1/2 h ~ 4m 2 Af| 7 - 1/2 /(*- 1 ) B 23fl , i^ 3 7 J ' fc V? ) P B • (3-24) 

Here,V^ p B = dkP B + T B D dk4> c 'p D , with T B D being the Levi-Civita connection of the metric 
Jab- To get (3.24), we have eliminated Tab using the 0> AB equation, 

Tab^^Vcd-O, (3.25) 

and in deriving the second term we have made use of (3.15). Note that (3.24) still contains 
the $ , but it is easy to see that it depends only on the 9 combinations of these that are 
determined in terms of other variables by ^ab = 0: Multiplying (3.9) by (B~ x ) k a (b a )c& CA , 
where B^ 1 is defined in (3.16), we obtain, 

(B- X ) k a (b a )c* CA (fAB + ^) = (S-^BDa^V* • (3-26) 

Hence the equation ^ ab = that determines only 9 components of § , contains this matrix 
in the same combination in which it appears in (3.24). Thus on gauge fixing the (j) by using 
coordinate transformations, solving for the p using the N and N' 1 equations and eliminating 
the 9 components of $ using ^>ab = 0, we obtain a constraint £2(7, 7r ) ~ 0. This has the 
desired form to eliminate the momentum canonically conjugate to the ghost field. Finally 
note that (3.25) along with (3.14) implies that Tab = ^PaPb- 
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